WAVE EQUATION SOLUTIONS AND PAIR PRODUCTION 
FOR ARBITRARY SPIN PARTICLES 
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Abstract 

We investigate the theory of particles with arbitrary spin and mag- 
netic moment in the Lorentz representation (0, s)©(s, 0) in an external 
constant and uniform electromagnetic field. We obtain the density ma- 
trix of free particles in pure spin states. The differential probability of 



CO ■ pair producing particles with arbitrary spin by an external constant 



and uniform electromagnetic field is found using the exact solutions. 
We calculate the imaginary and real parts of the Lagrangian in an 
electromagnetic field that takes into account the vacuum polarization. 
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^ 1 Introduction 

Interest in the theory of relativistic particles with arbitrary spins is increas- 
ing. One of the reasons is that SUSY models require superpartners, i.e., 



additional fields of particles with higher spins. In particular, it is important 
to take into account particle with spin 3/2-gravitino which appears in cos- 
mological models based on supergravity and in the theory of inflation of the 
universe [1]. 

Another factor is that some string models have similar features to models 
of relativistic spinning particles [2] . It is also interesting to have particles with 
arbitrary fractional spins [3] (see also refs. [4,5]). Such spinning particles 
in (2 + 1) dimensions, called anyons, were discovered and have anomalous 
statistics. 

There are many different relativistic wave equations which describe par- 
ticles with arbitrary spins [6-11]. The fields of free particles of a mass m and 
spin s in these formalisms realize definite representations of the Poincare 
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group. Some of these schemes are equivalent to each other. If the inter- 
action with external electromagnetic fields is introduced, approaches based 
on different representations of the Lorentz group become inequivalent. Most 
theories of particles in external electromagnetic fields have difficulties such 
as non-causal propagation [12], indefinite metrics in the second-quantized 
theory [13-15], etc. 

We proceed from the second-order relativistic wave equation for particles 
with arbitrary spin s and magnetic moment /i, on the basis of the Lorentz 
representation (0, s) © (s, 0), where © means the direct sum. In this scheme 
the Hermitianizing matrix r\ exists for arbitrary spin, as there are pairs of mu- 
tually conjugate representations (0, s) and (s, 0) [11]. Therefore we have here 
a Lagrangian formulation and the scalar product of wave functions, which 
is important for any quantum mechanical calculations. Such an approach 
avoids difficulties of other schemes and the corresponding wave function has 
the minimal number of components. This is a generalization of the Feynman- 
Gell-Mann equation [16] for particles with spin 1/2 to the case of higher spin 
particles which possess arbitrary magnetic moment. In the particular case 
of spin 1/2 we recover to the Dirac theory. If the normal magnetic mo- 
ment is considered, it leads to the approach of refs. [11]. Particles in this 
scheme propagate causally in external electromagnetic fields and this is a 
parity-symmetric theory with a Lagrangian formulation. Some equations on 
the basis of (0, s) © (s, 0) representations of the Lorentz group were stud- 
ied in refs. [17]. In refs. [11], the 2(6s + 1) -component, first-order matrix 
parity-invariant formulation of the equation for particles with arbitrary spin 
was considered. Then the author obtained the second-order equation for 
particles with "normal" magnetic moment. Starting with the second-order 
equation for particles which possess arbitrary magnetic moment, we go to the 
first-order wave equation with another representation of the Lorentz group; 
the algebraic properties of the matrices are the same as in the approach of 
refs. [11]. We find solutions of equations for free particles in the form of den- 
sity matrices (projection matrix-dyads) for pure spin states which are used 
for different electromagnetic evaluation of the Feynman diagrams. Such pro- 
jection matrix-dyads allow us to make covariant calculations without using 
matrices of the wave equation in a definite representation. 

The main purpose of this paper is to investigate solutions of wave equa- 
tions, pair production of arbitrary-spin particles by constant, uniform elec- 
tromagnetic fields, and vacuum polarization of higher spin particles. Consid- 
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ering one-particle theory, and obtaining the differential probability for pair 
production of particles with arbitrary spins, we avoid the Klein paradox [18, 
19]. As a particular case of spin- 1/2 and gyromagnetic ratio 2 particles, we 
arrive at the well-known result found by Schwinger [20], who predicted e + e~ 
pair production in a strong external electromagnetic field. This has now been 
realized by the development of power laser techniques. It should be noted 
that the pair production of particles by a gravitational field is also important 
for understanding the evolution of the universe near a singularity [21]. 

The probability of pair production of particles in external electromagnetic 
fields can be found on the basis of exact solutions of the wave equations [22] 
or the imaginary part of the Lagrangian [20]. We consider here both ap- 
proaches. Nonlinear corrections to the Maxwell Lagrangian of the constant 
uniform electromagnetic fields are determined from the polarization of the 
vacuum of arbitrary-spin particles. The problem of pair production of parti- 
cles with higher spins using the quasiclassical scheme (method of "imaginary 
time") was considered in ref. [23] and is in accord with our approach via 
relativistic wave equations. It should be noted that the quasiclassical ap- 
proximation has a restriction for the fields E, H <C m 2 /e when the process 
is exponentially suppressed. It means that the approach of ref. [23] is valid 
when the electromagnetic fields are not too strong, i.e., less than the critical 
value m 2 /e. But it is known that pairs of particles are created rapidly at this 
critical value of the fields. In our consideration, there are no such restrictions. 
The problem of the pair production of vector particles with gyromagnetic ra- 
tio 2 was investigated in ref. [24]. In refs. [25, 26], the imaginary part of the 
effective Lagrangian which defines the probability of e + e~-production was 
found by taking into account the anomalous magnetic moment. 

We use system of units h — c — 1, a — e 2 /4n = 1/137, e > 0. In Sec- 
tion 2, proceeding from the second-order equation for arbitrary-spin particles 
with anomalous magnetic moment, we go to the first-order formulation of the 
theory. All independent solutions of the equation for free particles are found 
in the form of matrix-dyads (density matrices). In Section 3, we consider 
the important case of spin-1 particles, which is different from the Proca or 
Petiau-Duffin-Kemmer theories. Here, the 2 (6s + 1) -component first-order 
equation is constructed. Section 4 investigates exact solutions of arbitrary- 
spin particle equations in constant, uniform electromagnetic fields. On the 
basis of exact solutions, we find the differential probability of pair production 
of particles with arbitrary spin and anomalous magnetic moment. The imag- 
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inary part of the effective Lagrangian for electromagnetic fields is calculated. 
In Section 5, using the Schwinger method, we find the nonlinear corrections 
to the Lagrangian of a constant, uniform electromagnetic field caused by the 
vacuum polarization of particles with arbitrary spin and magnetic moment. 
Section 6 discusses of results. 

2 Wave Equation and Density Matrix 

We proceed here from the theory based on the (s, 0) © (0, s) Lorentz repre- 
sentation for massive particles. The wave function of the (s, 0) © (0, s) rep- 
resentation has 2(2s + 1) components. For spin 1/2, we arrive at well-known 
Dirac bispinors. For spin 1, however, there is doubling of the component 
compared with the Proca theory [27] because the vector particles have three 
spin states with the projections s z — ±1, 0. 

We postulate the following two (for e — ±1) wave equations for arbitrary- 
spin particles in external electromagnetic fields: 

- m 2 - ||i^Eg) * e (x) = 0, (1) 

where s is the spin of particles, = — ieA^; F^ v = d^A v — d v A^ is the 
strength tensor of external electromagnetic fields, e — ±1; and Yk£) 
are the generators of the Lorenz group which correspond to the (s, 0) and 
(0, s) representations. Two equations (1) (for e = ±1) describe particles 
which possess the magnetic moment \i = eq/(2m) and gyromagnetic ratio 
g = q/s. At q = 1, we have the "normal" magnetic moment \i = e/(2m) 
and g = 1/s. The generators SjfJ are connected with the spin matrices Sk 

by the relationships = tijkSk and = —ieSk, where the parameter e 
corresponds to the Lorentz group representations (s, 0) for e = 1 and (0, s) 
for e = — 1. As usual, relations 

[Si, Sj] = ie ijk S k , (S 1 ) 2 + (S 2 ) 2 + (S 3 ) 2 = s(s + 1) (2) 

are valid, where i,j, k = 1,2, 3; is the Levi-Civita symbol (ei23 = 1). 

At q — 1 equations (1) where considered in refs. [11]. The theory of 
arbitrary spin particles based on Eq. (1) is causal in the presence of ex- 
ternal electromagnetic fields. It is seen from the method of ref. [12] that 
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the equations remain hyperbolic and the characteristic surfaces are lightlike. 
Equation (1) is invariant under the parity operation. Indeed, at the parity 
inversion e — * —e, the representation (s, 0) transforms into (0, s). 

Now we will consider the problem of formulating a first-order relativistic 
wave equation from the second-order equation (1). This is convenient for 
some quantum electrodynamic calculations with polarized particles of arbi- 
trary spins. 

Let us introduce the matrix e A)B with dimension nxn; its elements consist 
of zeroes and only one element is unity, where row A and column B cross. 
So the matrix elements and multiplication of this matrices are 

(c A ' B ) CD = SacSbd, e A > B e c > D = 6 BC e A > D , (3) 

where indexes A, B,C,D — 1, 2...n. 

Six generators Ej^ (or Sj^) occurring in Eq. (1) have the dimension 
2s + 1. Therefore the wave function ^+(x) (or of Eq. (1) possesses 

2s + 1 components. Now we can introduce the 5(2s + l)-component wave 
function 

•'«-(-±5&w) (4) 

so that = {^a(x)}, A = 0,/i; ^ = ^ = -{l/m)D^ + (x) 

and ^+(x) realizes the Lorentz representation (s, 0). It is not difficult to 
check that Eq. (1) for e = +1 can be represented as the first-order equation 



/^D ll + m)V 1 (x) = 0, (5) 
where 5(2s + 1) x 5(2s + 1)- matrices 

(3^ = (e°>» + e* ) ® I 2s+1 - £e * ® Ej+), (6) 

are introduced, and <S> is the direct product of matrices, I2S+1 is a unit matrix 
of the dimension 2s + 1, and in (6), we imply the summation on index v. 
It should be noted that in refs. [11], the (6s + 1) -component wave func- 
tion was introduced for the case of "normal" magnetic moment of particles. 
Here the higher dimension 5 (2s + 1) of the wave function is considered as 
compared with the (6s + 1) -component function investigated in refs. [11]. 
The difference is that we do not introduce the vector potentials for arbitrary 
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spin fields and deal only with the strength field tensor ^+(x) and its deriva- 
tives D^ + (x). The particular case of a spin-1 field will be considered in 
Section 3 on the basis of (6s + 1 ^representation of the wave function by in- 
troducing vector potentials. Using properties (3), it is easy to check that the 
five-dimension matrices P^ DK = e 0,fl + obey the Petiau-Duffin-Kemmer 
algebra [28-30], 

nPDK nPDK nPDK , nPDK nPDK nPDK r nPDK , r nPDK / 7 \ 
P M Pv Pa +Pa Pv Pn = <VA* + <W/^ . (7) 

The wave function transforms on the [(0,0) © (1/2,1/2)] ® (s, 0) 

representation of the Lorentz group [31,32]. For the case e = — 1, we have 
the analogous equation 



where 



(/3H^ + m)* 2 (a;) = 0, (8) 



= ( -i Wx) ) (9) 

and * 2 (x) = {^b(^)}, S = 0,71; %(ar) = ^ = -^L^M^), where 

transforms as (0, s) representation of the Lorentz group, and ^2 0*0 
realizes the [(0,0) © (1/2, 1/2)] © (0, s) representation. The two equations 
(5) and (8) can be combined into one first-order equation 

(P II D II + m)^(x) = (10) 
with the matrices and wave function 

P, = P^®P!r\ ^ = (tlU)- < n ) 

Using the properties of elements of the entire algebra (3), it is not difficult 
to verify that the matrices (3^ (the same for (3^ and Pj~^) obey the algebra 

PpPvPa + PvPoP» + PaP^Pv + PvP^Pa + PpPaPv + PaPvP^ = 

= 2 (<WA, + S^p v + S^p a ) . (12) 
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In other refs. [11], the 2(6s+l)-dimensional representation of the SL(2, C) 
group was considered for particles of arbitrary spin with algebra (12). Here 
we study another representation of the SL(2, C) group, and the matrices (3^ 
of (11) are 10 (2s + 1) x 10 (2s + 1) dimensional. Different representations of 
this algebra were studied in refs. [33] and [34]. 

Let us consider the problem of finding the solutions to Eq. (10) for defi- 
nite momentum and spin projection. It is convenient to find these solutions 
in the form of projection matrix-dyads (density matrices). All electrody- 
namic calculations of Feynman diagrams with arbitrary-spin particles can be 
done using these matrices. As particles in initial and final states are free 
particles, we can put the parameter q = in (1), (10). This corresponds to 
the case when external electromagnetic fields are absent. Then the matrices 

transform to : 

$1 = [(e^ + e*°) <g> J 2s+1 ] © \(e^ + <g> I 2s+1 )] , (13) 

which obey the Petiau-Duflin-Kemmer algebra (7). The projection operators 
extracting states with definite 4-momentum p^ for particle and antiparticle 
are given by 

M± = !£ (|±™) (14) 

where p = (we use the metric such that p 2 = p 2 +p\ = p 2 — Pq = —m 2 ). 
Plus and minus signs in (14) correspond to the particle and antiparticle, 
respectively. Matrices A± have the usual projection operator property 



Ml = M±. (15) 

Equation (15) is verified by the relation p 3 = p 2 p, which follows from the 
Petiau-Duffin-Kemmer algebra (7). To find the spin projection operators, we 
need the generators of the Lorentz group in the representation of the wave 
function fy(x) which occurs in Eq. (10). From the structure of the functions 
^i(x) and ^2(x) in (4) and (9), we define the generators of the Lorentz group 
in our 10(2s + l)-dimensional representation 

/ = /(+) ffi /(-) 
j(+) = _ / 2s+1 + ih £(+), (16) 
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jp = (?* - ^) ® w + ih ® , 

where J5 is the 5-dimension unit matrix. Using properties (3), we get the 
commutation relations for J^, and 

[Jfiu, Ja/3] = 5 va J[t/3 + S^Jua — 5^ a J v [3 — 5 u /3.J fMa , (17) 
[Pfji, Jap] = — S^pPa- (18) 

Relationship (17) is a well-known commutation relation for generators 
of the Lorentz group [31,32]. Equation (10) is a form invariant under the 
Lorentz transformations because relation (18) is valid. To guarantee the 
existence of a relativistically invariant bilinear form 

*(x)*(a:) = ^ + (x)r]^(x), (19) 

where ^/ + (x) is the Hermitian-conjugate wave function, we should construct 
a Hermitianizing matrix r] with the properties [33,35]: 



riPi = -PiV, Vfa = far) (i = 1,2,3). (20) 

Such a matrix exists because here there are two (0, s) and (s, 0) mutually 
conjugate representations; 77 is given by 

v = ( e °> s + £ ^ a - e 4 ' 3 - e 3 ' 4 - e ' 5 - e°>°) <g> I 2s+1 , (21) 

where the summation on the index a = 1, 2, 3 is implied. Now we introduce 
the operator of the spin projection on the direction of the momentum p: 

^ % 

S P = ~7T\ \ e abcPaJbc = («p + <T p ) © (~K p + W p ) , (22) 

where 

£ 2 — _ 

/%> = - | -^abcPa£ b,C ® hs+l, Kp = ~] ^abcVa^ ® hs+l, 
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o p = Op = h ® y~~~y ( 23 ) 

and | p \ = \Jp{ + v\ + V%- It is eas Y to check that the commutation relation 
holds: 



= 0. 



The matrices k p , k p obey the simple equations 

k p = K pi K p = K p- (24) 

Taking into account Eqs. (2), we derive relations for the matrices a p (23): 

^(Tp — • ■ • (a 2 — s 2 ) = for half-integer spins, 

°"p (o 2 — l) • • • ((Tp — s 2 ) = for integer spins. (25) 

Relations (25) allow us to construct projection operators which extract 
the pure spin states. Using the relationship 

S P P = (o p © Op) p, 

we can consider projection matrices on the basis of equations (25). The 
common technique of the construction of such operators is described in refs. 
[35]. Let us consider the equation for the auxiliary spin operators a p , a p for 
the spin projection s k : 

a p ^ k = s k ^ k . (26) 

The solution to Eq. (26) can be found using relationships (25), which can 
be rewritten as 

(v P - s k ) P k (a p ) = 0, (27) 
where the polynomials P k (cr p ) are given by 

Pk(cr p ) = (a p - • • (a p + s k ) ■ ■ ■ (<7 2 - s 2 ) for half-integer spins, 
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Pk(o-p) = (Jp (a* - l) • • • (cr p + s k ) ■ ■ ■ (crp - s 2 ) for integer spins. (28) 

Every column of the polynomial P k {(J p ) can be considered as an eigenvec- 
tor ^/ k of Eq. (26) with the eigenvalue s k . As s k is one multiple root of Eqs. 
(25), all columns of the matrix P k {a p ) are linearly independent solutions of 
Eq. (26) [35]. It can be verified using definitions (28) that the matrix 

Pk{Sk) 

is the projection operator with the relation 

Ql = Q k . (30) 

Equation (30) tells that the matrix Q k can be transformed into diagonal 
form, with the diagonal containing only ones and zeroes. So the Q k acting 
on the wave function \1/ will retain components which correspond to the spin 
projection s k . 

We have mentioned that this theory of arbitrary-spin particles involves the 
doubling of the spin states of particles because there are two representations 
(s, 0) and (0, s) of the Lorentz group. To separate these representations, 
which are connected by the parity transformations, we use the parity operator 

K = (e^ + + + e°>°) <g> I 2s+1 (31) 



with the summation on index /j, — 1,2, 3, 4. The 10(2s+l)-dimensional matrix 
K has the property K 2 = iio(2s+i)- The projection operator extracting states 
with the definite parity is given by 

\ £ = l -(K + s), (32) 
where e — ±1. This matrix possesses the required relationship 

A 2 = A e . (33) 

It should be noted that the matrix K (31) plays the role analogous to the 
7 5 matrix in the Dirac theory of particles with the spin 1/2. It is easily verified 
that the operators p, S p , K commute with each other, and, as a consequence, 
they have common eigenvectors. The projection operator extracting pure 
states with definite 4-momentum projection, spin, and parity is given by 
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Ii ±mAe = M±A £ (Q k © Q k ) (34) 

with matrices (14), (29) and (32). The H± m ,k,e is the density matrix for pure 
states. It is easy to consider impure (mixed) states by summation of (34) 
over definite quantum numbers Sk, e. The projection operator for pure states 
can be represented as matrix-dyad [36] 

n ±m , fc , £ = n m , k • % m ,k, (35) 

where \&± mjfc = {^± m ,k} + Vi an d ^± m ,k * s the solution to equations 



Kn m , k = e*W (36) 

Expression (35) is convenient for calculations for different quantum elec- 
trodynamics processes with polarized particles of arbitrary spins. 



3 Wave Equation for Spin-1 Particles 

As a particular case, we consider the wave equation for particles with spin 
I. To compare our approach with refs. [11], the 2 (6s + 1) representation for 
the wave function will be studied. Equation (1) in the case of spin-1 particles 
becomes 

(p\ - m 2 ) ^^(x) - ieq5[ llu ][ a fs]F aa ^ a f 3 (x) = 0, (37) 

where wave function ^ ^ u {x) is the antisymmetric tensor of the second rank 
^ \tu( x ) = —^vn(x), which has six independent components, and 5[^„][ aj g] = 
— bupbva- For comparison with refs. [11], the "normal" magnetic 
moment \i = e/2m will be considered here. Then the parameter q — 1, and 
gyromagnetic ratio g — 1. Let us introduce two 4- vector potentials ip^x) 
and ipfj,(x) in accordance with the relationship [36] 

V^ix) = D^ v (x) - Dyip^x) - e^apD^^x), (38) 

where e^ va p is an antisymmetric tensor Levi-Civita, €1234 = —%■ This is a 
more general representation of the antisymmetric tensor \1/ ^ix) via 4- vector 
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ipn{x) and 4-pseudovector ip^{x). It is not difficult to verify that Eq. (37) 
(at q — 1) is a consequence of equations 

D^^x) + m 2 VvO*0 = °> ( 39 ) 

D u ^^(x)+m 2 4,^x) = 0, (40) 

where ^^(x) = ^^^^^(x) is the dual tensor. So, the system of the first- 
order equations (38)-(40) is equivalent to the second-order equation (37). 
A doubling of spin states of particles is obvious here. It follows from Eqs. 
(39) and (40) for free particles (when external electromagnetic fields are 
absent and = — ieA^ — > <9 M ) that the Lorentz equations d^ip^x) = 0, 
dfj,ipfj,(x) = hold. Therefore the vector potential ift^x) describes spin-1 
states (the state with spin is not present due to the Lorentz equation) and 
pseudovector ip^x) also describes spin-1 states. As a result, there are six 
spin states and the system of equations (38)- (40) is not equivalent to the 
Proca theory [27]. If we set ipn(x) = in (38)-(40), we will arrive at the 
Proca equations. 

Equations (38)-(40) are equivalent to the following system: 

D v M„v + m 2 M^ = 0, 

M M „ = D^M U - D U M^ - e^pDaMp, (41) 
with the self-dual tensor M M „ = iM^ v , and 

D U N^ + m 2 N^ = 0, 

Nfu, = D^N U - D U N^ + e^apDnNp, (42) 
with the ant i- self- dual tensor = —iN^ v , where 

= X - (%{x) - i^ix)) , = X - + M^(x)) , 

= - (M x ) + iM*)) > N »v = g (V(^) - ^^(x)) . 

Adding and subtracting Eqs. (41), (42), we get Eqs. (38)-(40). The 
self-dual tensor M^ v which obeys Eqs. (41) is transformed under (1,0) rep- 
resentation of the Lorentz group and has three independent components. 
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Equations (41) are not invariant to the parity transformation and there is 
no Lagrangian formulation of them. This also applies to Eqs. (42) for the 
ant i- self- dual tensor N^ u , which transforms under (0, 1) representation of the 
Lorentz group. But if we consider the hole system of equations (41), (42) 
[which is equivalent to Eqs. (38)-(40)] on the basis of (1,0) © (0,1) repre- 
sentation of the Lorentz group, we will have P-invariant theory within the 
Lagrangian formulation. 

Consider now the first-order formulation of Eqs. (38)-(40) in matrix form 
(10). Introducing a 14-dimensional wave function 

(if) (x) \ 
C(z) , (43) 

where index A — //, Jl so that * M (x) = W( x )> ^^(x) = ^^(x), 

ty~(x) = Vv(x), we have that Eqs. (38)-(40) take the form (10) with the 
matrices 

h = PP+PF>, (44) 
i g(i) =e W +£ [A(i] l A ) ( 45 ) 

= (e~ XM + e [a ^) , (46) 

where [fiu] means the antisymmetric combination of indexes fi and v and 
corresponds to the six-dimensional subspace and there is a summation on 
the repeating indexes A in (45), (46). Here the representation with the di- 
mension 2(6s + 1) = 14 is valid. The (6s + 1) representation is built from 
and self-dual tensor M^, or from and the anti-self-dual tensor N^. The 
resulting 2 (6s + 1) representation (43) is the direct sum of the above repre- 
sentations. The fact that matrices (44) obey the algebra (12) is confirmed 
by the equalities (3). This representation is different from the 10(2s + 1)- 
dimensional representation (11). For spin-1 particles, the wave function (11) 
has 10 (2s + 1) = 30 components and is given by 

The approach considered in Section 2 allows us to apply the unique tech- 
nique for arbitrary-spin particles which possess arbitrary magnetic moment 
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without introducing "potentials" , but the cost of this is the high dimension 
10 (2s + 1). 

It should be noted that 10- dimensional matrices f3^ in (45) obey the 
Petiau-Duffin-Kemmer algebra (7), and that they act in 10-dimensional sub- 
space of wave functions 

*"*<*> = ( tw ) • 

where we imply that ^ ^{x) = D fM ip u (x) — D^ift^x), so that 10-dimensional 
Petiau-Duffin-Kemmer equation is given by 

(pWD fl + m)v PDK (x)=0. (47) 

This corresponds to spin-1 particles without the doubling of spin states. 
The generators of the Lorentz group in our 14-dimensional representation 
(43) are given by 

= £ ^ - £ ^ + £ M,[*"] _ £ M.M + ~w _ £ v,t. ( 48 ) 

they obey the required commutation relations (17), (18). The projection 
operators extracting states with definite 4-momentum have virtually the 
same form (14), but with matrices (44). The operator for the spin projection 

S P = _ 771 i e abcPaJbc, (49) 

2 | p | 

with generators (48) satisfies the equation 

(Sp - l) = 0. (50) 

In accordance with this procedure [see (27)- (30)], we find the projection 
operators corresponding to the definite spin projections 

Q± = \s p (S p ± l) , (51) 

Qo = 1 - S 2 p . (52) 
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Operators Q±, Qo extract the spin projections Sk = ±1 and Sk = 0, 
respectively. Equations (38)- (40) are invariant under the dual 5*0(2) trans- 
formations 

ip'^x) = cosa + ^^{ x ) sin a, 

ip'^x) = —ipn(x) sin a + ip^(x) cos a. (53) 

This group of symmetry is connected with the doubling of spin states of 
particles. To remove this degeneracy, we use the projection generator (32) 
with the matrix 

K = + + ^ Q/3£ [H,M (54) 

which has the property K 2 = I14, where 1 14 is the unit matrix in the 14- 
dimensional space. Now we construct the density matrices of spin-1 particles 
for pure states as products of the projection operators: 

(S P + s k ) ip {ip ±m)(K + e), (55) 

n ±mAe = (l - S 2 p ) ip{ip±m) (K + e), (56) 

where p = and the spin projections on the direction of the momentum 

p are Sk = ±1 in (55) and Sk = in (56), and e = ±1. Matrices (55), (56) 
can be represented in the form of matrix-dyads (35) with the Hermitianizing 
matrix 

^ _ £ m,m _ £ A,4 _|_ £ [mA],[mA] _ \ £ [mn] ,[mn] _ £ m,m _|_ £ 4,4 

2 

where m, n = 1,2,3, which obeys the general equations (20) with matrices 
(44). The Lagrangian of free particles (where the vector potential of the 
electromagnetic field = 0) takes the standard form of 

C = -*(x) (f3^ + m)V(x), 

where ^/(x) = ^ + (x)i], and ty + (x) is the Hermitian-conjugate wave function. 
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4 Pair Production by External Electromag- 
netic Fields 

To calculate the probability of pair production of arbitrary-spin particles, we 
follow the Nikishov method [22]. Thus, exact solutions to Eq. (1) should 
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be found for external, constant, uniform electromagnetic fields. Using the 
properties of generators we find the relationships 

2^v F ^ — SiXi, -T^F^ = SiX*, (58) 

where Xi = H \ + iEi, X* = Hi — iEf, Ei, Hi are the electric and magnetic 
fields, respectively, and the spin matrices Si obey Eqs. (2). In the diagonal 
representation, the equations for the eigenvalues are 



SiXM+\x) = aXV { +\x), SiX*¥°\x) = <tX*¥?\x), (59) 

where X = VX?, X = H+iE, and the spin projection a is 



±s, ±(s — 1), • • -0 for integer spins, , , 

±s, ±(s — 1), • • • ± \ for half-integer spins. 

Taking into account (58) and (59), we rewrite Eq. (1) (for e = ±1) as 



(Dl-m 2 -aaX)¥°>(x)=0, (L> 2 - m 2 - aaX*)¥y(x) = 0, (61) 

where a = eq/s. These equations are like the Klein-Gordon equation for 
scalar particles, but with complex "effective" masses m 2 ^ = m 2 + aaX and 

(mljrj^j = m 2 + aaX*. It is sufficient to consider only one of Eqs. (61). Let 
us consider the solution of the equation 

{Dl - m 2 eff )^\x) = 0, (*<">(*) = *%\x)) (62) 

in the presence of the external, constant, uniform electromagnetic fields. 
The general case is when two Lorentz invariants of the electromagnetic fields 
F = \ F lv 7^ and £ = jF^F^ ^ (F^ = |e^ a/3 F a/3 ; e^ af3 is the 
antisymmetric Levi-Civita tensor). Then there is a coordinate system in 
which the electric E and magnetic H fields are parallel, i.e., E || H. In this 
case, the 4-vector potential is given by 

— (0, X\H, —xqE, 0) , (63) 
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so that 3-vectors E = nE and H = nH are directed along the axes, where 
n = (0, 0, 1) is a unit vector. The four solutions of Eq. (62) for the potential 
(63) with different asymptotic forms are given in refs. [22, 37] (see also refs. 
[38]) 

= iVexp [i(p 2 x 2 + p 3 x 3 ) - ^ | H n ( V )^\r), (64) 
where N is the normalization constant, H n {rj) is the Hermite polynomial, 



eHxi + p 2 ik 2 1 

r/ = 

and 



V V^H ' V 2eE 2' 



T=^(x + ^) 



+ ^\t) = D v [-(1 - i)r], ~^\r) = D u [(l - i)r], 

+^W(r) = D u *[(l + i) T ], _^(r) = D v * [—(1 + i)r]. (65) 

Here D u (x) is the Weber- Hermite function (the parabolic-cylinder func- 
tion). The probability for pair production of particles with arbitrary spins by 
constant electromagnetic fields can be obtained through the asymptotic form 
of the solutions (65) when the time xo — > ±oo. At x — > ±oo, the functions 
+*0( ct )(t) have positive frequency and Zip^(r) have negative frequency. The 
constant k 2 which enters the index v of the parabolic-cylinder functions (65) 
is given by [38] 

k 2 = m 2 eff + eH(2n + l), (66) 

where n — I + r, I is the orbital quantum number, r is the radial quantum 
number, and n — 0, 1,2, ... is the principal quantum number. It should be 
noted that for scalar particles, we have the equation k 2 = p$—p\, where po is 
the energy and p 3 is the third projection of the momentum of a scalar particle. 
In our case of arbitrary-spin particles, the parameter m 2 ^ is a complex value. 
Nevertheless all physical quantities in this case are real values. Solutions (64), 
(65) are characterized by three conserved numbers: k 2 and the momentum 
projections p 2 , Pz- As shown in refs. [22], the functions (64) are connected 
by the relations 

= c lnr7 + ^l(x) + c 2na -^ n (x), 
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where coefficients Ci n<7 , c 2no - are given by 



(67) 



C2r, 



exp 



11 ,x -\ 



_ m 2 eff + eH(2n + l) 



|2 _ I |2_ 



for bosons, 



+ c 2r , 



for fermions. 



(68) 



The values Ci na , c 2n o- are connected with the probability of pair production 
of arbitrary-spin particles in the state with the quantum number n and the 
spin projection a. The absolute probability for the production of a pair in 
the state with quantum number n, momentum p and the spin projection a 
throughout all space and during all time is 



C2na | = exp 



-7T 



771 

eE 



qaH H 



sE 



E 



(2n 



(69) 



The value (69) is also the probability of the annihilation of a pair with 
quantum numbers n, p, a with the energy transfer to the external electro- 
magnetic field. It is seen from (69) that for H 3> E, the pair of particles 
are mainly created by the external fields in the state with n = 0, a = —s. 
This is the state with the smallest energy. So at H ^> E there is a produc- 
tion of polarized beams of particles and antiparticles with the spin projection 
a = — s ( s is the spin of particles). The average number of pairs of particles 
produced from a vacuum is 



C2r 



dp 2 dp 3 



L 2 



(70) 



because {2ti)~ 2 dp2dp 3 L 2 is the density of final states, where L is the cutoff 
along the coordinates, i.e., L 3 is the normalization volume. The variables 77, 
r define the region where the process occurs, which is described by solutions 
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(64) with the coordinates of the center of this region xq = —ps/eE, x\ = 
—p2/eH. Therefore instead of the integration over p 2 and ps in (70), it is 
possible to use the substitutions [22] 

J dp 2 -> eHL, J dp 3 -> eET (71) 

with T being the time of observation. 

Evaluating the sum in (70) over the principal quantum number n, with 
the help of (69), (71), we obtain the probability of pair production per unit 
volume and per unit time 

. , N e 2 EHexp\-nm 2 /(eE)}^ . 

where 6 = 7iqH/(sE), V = L 3 . Now we evaluate the sum over the spin 
projection a in (72) for integer and half-integer spins: 
1. Integer spins: 

-b(s+l) _ 1 

J2exp{-ab) =S! + S 2 , S 1 = e° + e~ b + ... + e~ sfe = — - b 



e~ b - 1 



9 „» + ^4. + ^ 6 ( e - 1 ) c | c cosh( g fe) - cosh [(s + 1)6] 
S2-e+e +... + e - , S, + S 2 x _ COBhb > 

(73) 

2. Half-integer spins: 

p -b(s+l) _ -6/2 

E exp(-<7&) = ^1 + S 2 , S[ = e^ 2 + e-^ + ... + e- sb = , 

S' = e b ' 2 +e 3b l 2 + +e sb = — ~ — S'+S' = cosh (*ft) ~ cosh [(s + 1)6] 
2 "" e b -l ' 1 2 1- cosh b 

(74) 

So the final expressions for integer and half-integer spins, (73) and (74), 
are the same. Using the relationship 

cosh(sfr) - cosh Us + 1)6] , . . , . ,6 sinh \b(s + 1/2)] 

^ l i = cosh (sb) + sinh(sfe) coth - = 1 ; - , / n 

1- cosh b K J K J 2 sinh(6/2) 

(75) 
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and equations (72), (73) we arrive at the pair-production probability 



HE H) = — = e2EH 6XP [~ 7rm2 /( eE )] sinh K 2s + ^H/(2sE)} 
1 ' ' VT ' 8tt 2 sinh (irH/E) sinh [qirH/(2sE)} ' 1 ' 

Expression (76) coincides with that derived in refs. [23] using the quasi- 
classical approach. So I(E, H) is the intensity of the creation of pairs of 
arbitrary-spin particles which possess the magnetic moment /i = eq/(2m) 
and gyromagnetic ratio g — q/s. 

In ref. [23], there is a discussion of physical consequences which follow 
from Eq. (76). In particular, there is a pair production in a purely magnetic 
field if q — gs > 1 [23]. It is interesting that the exact formula derived here 
from quantum field theory (which is valid for arbitrary fields E, H) coincides 
with the asymptotic expression obtained in ref. [23] for E,H <C m 2 /e. 

To obtain the imaginary part of the density of the Lagrangian, we use 
the relationship [22] 

VTlmC = \ [ ]T In | c lna | 2 d P2 d P3J f^. (77) 

^ J n,cr K^) 

With the help of (68), (71), we arrive at (see also ref. [23]) 



T „ e 2 EH ~ f3 n ( Tim 2 n\ sinh [n(2s + l)qnH/(2sE)] 



where 



16-7T 2 ~^ n \ eE J sinh (nirH/E) sinh [nqnH/ (2sE)} ' 

(78) 

(— l) n_1 for bosons, 

1 for fermions. 



The different expressions for bosons and fermions occur due to different 
statistics and relations (68). The first term (n — 1) in (78) coincides with 
the probability of the pair production per unit volume per unit time divided 
by 2 [22] (see discussion in ref. [23]). 
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5 Vacuum Polarization of Arbitrary-Spin Par- 
ticles 



Now we calculate the nonlinear corrections to the Lagrangian of a constant, 
uniform electromagnetic field interacting with a vacuum of arbitrary-spin 
particles with gyromagnetic ratio g. For the case of spins 0, 1/2 and 1 (for g = 
2) this problem was solved in refs. [39,40,20,24]. The nonlinear corrections to 
the Lagrangian of the electromagnetic field describe the effect of scattering of 
light by light. We consider one- loop corrections to the Maxwell Lagrangian 
corresponding to arbitrary-spin particles, and, to take into account vacuum 
polarization, it is convenient to adapt the Schwinger method [20]. Applying 
this approach to the arbitrary spin particles described by Eq. (1), we arrive 
at the effective Lagrangian of constant, uniform electromagnetic fields, 



, v;i 2 J o drr 3 exp (-m 2 r - Z(r)) tr exp ^E^f^rj , (79) 

where e = 1 for bosons, e = — 1 for fermions, 

5V = E« © £<->, J(r) = itrln [(eFr)" 1 sin(eFr)] (80) 

and F^y is a constant tensor of the electromagnetic field. The formal dif- 
ference of (79) from the case of spin-1/2 particles is in the substitution 
°> ~~ *• (<?/ s )£W> where = (i/2) [7^, 7„] , 7^ being the Dirac matrices. 
The problem is to calculate the trace of the matrices occurring in the expo- 
nential factor in (79). Using relations (58)-(60) and (73)-(75), we find 

trexp (^-H^F^T^j = 2Re cosh(egXr) + sinh(egXr) coth — ^ • 

(81) 

Inserting (81) into (79) and adding the constant which is necessary to 
cancel £1 when the electromagnetic fields are turned off [20], we arrive at 

£1 = — - / c?TT~ 3 exp (— m 2 T) X 
8n z Jo V ' 

(erfg 



x 



2 Re [cosh(egXr) + sinh(egXr) coth (eqXr/(2s))} 2s + 1 



2Imcosh(eXr) 



(82) 
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where Q = EH. With q — 1 and s = 1/2, this Lagrangian (82) coincides 
with that of Schwinger [20]. Expression (82) is the correction to the Maxwell 
Lagrangian that takes into account the vacuum polarization of arbitrary- spin 
particles which possess the magnetic moment fi = eqj (2m) and gyromagnetic 
ratio g = q/s. Adding (82) to the Lagrangian of the free electromagnetic fields 

Co = -T =\ (E 2 - H 2 ) 

and introducing the divergent constant for weak fields, we get the expres- 
sion for the total Maxwell Lagrangian 

C M = C + Ci = -ZT + ^— ^ J drr~ 3 exp [-m 2 Tj x 



A2 ^ Re [cosh(egXr) + sinh(egXr) coth (egXr/(2s))] 2s + 1 Q( 2 

' ieT) y 2Imcosh(eXr) 2 * 

(83) 

where 



te 2 /3 r . / 2 \ , 3 (2s + 1) [s (<J 2 - 1) + <J 2 ] 



24s 

(84) 

Schwinger's procedure is used to renormalize the electromagnetic field 
T — > ZjF and the charge e — > Z _1 / 2 e. After expanding (84) in the small 
electric and magnetic H fields, we arrive at the Lagrangian of a constant, 
uniform electromagnetic field (in rational units) 



C M = i ( E 2 - H 2 )-|^ [(E 2 - H 2 ) 2 (15,5 - 7 ) + (EH) 2 (47 + *±±)j +. 

(85) 

where a = e 2 / (47r) and 

= (2s+l)[g 4 (^ + l)(3 g 2 + 3s-l)-3s 3 ] _ 
16s 3 

It is easy to verify that for the particular case of s = 1/2, q = 1 (which 
corresponds to the Dirac theory), Eq.(85) coincides with the well-known 
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Schwinger Lagrangian [20]. With s = 1 and q = 2 expression (85) is dif- 
ferent from one obtained in ref. [24]. This is because our theory of particles 
with s — 1 is not equivalent to the Proca theory (see Section 3); there is 
the doubling of spin states here. The effective Lagrangian (85) is of the 
Heisenberg-Euler type, which has been found for the case of the polarized 
vacuum of particles with arbitrary spins and magnetic moment. Here we took 
into account virtual arbitrary-spin particles, but not virtual photons. This 
is because at small energies of the external fields, the radiative corrections 
are small quantities. It is not difficult to find the asymptotic form of (83) 
for supercritical fields eE/m 2 — > oo and eH/m 2 — > oo. It should be noted, 
however, that for strong electromagnetic fields, the anomalous magnetic mo- 
ment of electrons depends on the external field [41,42] and hence there is 
a similar dependence for arbitrary-spin particles. Therefore, to obtain the 
correct limit, it is necessary to take into account this dependence [26]. 

It is possible to obtain the imaginary part of the Lagrangian (78) from 
(83) using the residue theorem, taking into account of the poles of expression 
(83) and passing above them [20]. 

6 Discussion of Results 

The theory of particles with arbitrary spins and magnetic moment based on 
Eq. (1) and the corresponding Lagrangian allow us to find density matrices 
(34), (35), (55), and (56), the pair-production probability (76), and the effec- 
tive Lagrangian for electromagnetic fields (85), taking into account the po- 
larization of the vacuum. It is convenient to use matrix-dyads (34), (35), (55), 
and (56) for different electrodynamic calculations in the presence of particles 
with arbitrary spins. The exact formula for the intensity of pair production 
of arbitrary-spin particles coincides with the expression obtained in ref. [23] 
using the quasi-classical method of "imaginary time" , which is valid only for 
E, H <C m 2 /e, i.e., for weak fields. Hence it follows that the analysis in ref. 
[23] is valid for arbitrary electromagnetic fields and that it is grounded in rel- 
ativistic quantum field theory. In particular, there is a pair production by a 
purely magnetic field if gs > 1 [23], and in the presence of the magnetic field 
the probability decreases for scalar particles and increases for higher spin 
particles. As all divergences and the renormalizability are contained in Re£ 
(83), but not in Im£, the pair- production probability does not depend on the 
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renormalization scheme. The vacuum polarization corrections for scalar and 
spinor (with g — 2) particles are reliable because their theories are renormal- 
izable. The general formula (83) is of interest for the further development 
of the field theory of particles with higher spins (see discussion in refs. [23] 
and [24]). Expression (83) is a reasonable result for arbitrary values of s 
and q because for the particular case of scalar and spinor particles, we arrive 
at known results. Thus, we have here a reasonable and noncontradictory 
description of the nonlinear effects that arise in this interaction. 

References 

[1] G. Felder, L.Kofman, A. Linde, Phys. Rev. D60, 103505 (1999). 

[2] I. L. Buchbinder and SH. M. Shvartsman, Int. J. Mod. Phys. A8, 683 
(1993). 

[3] F. Wilczek, Fractional Statistics and Anyon Superconductivity (World 
Scientific, 1990). 

[4] J. L. Cortes and M. S. Plyushchay, Int. J. Mod. Phys. All, 3331 (1996). 

[5] S. I. Kruglov, Mod. Phys. Lett., A12, 1699 (1997); Europhys. Lett. 44, 
361 (1998); 

S. I. Kruglov and M. N. Sergeenko, Mod. Phys. Lett. A12, 2475 (1997). 

[6] H. J. Bhabha, Rev. Mod. Phys. 17, 200 (1945). 

[7] M. Fierz and W. Pauli, Proc. Roy. Soc. (London) A173, 211 (1939). 

[8] W. Rarita, J. Schwinger, Phys. Rev. 60, 61 (1941). 

[9] V. Bargman, E. Wigner, Proc. Am. Acad. Sci. 34, 211 (1948). 

[10] M. A. K. Khalil, Progr. Theor. Phys. 58, 1538 (1977). 

[11] W. J. Hurley. Phys. Rev. D4 3605 (1971); D10 1185 (1974); Phys. Rev. 
Lett. 29, 1475 (1972). 

[12] G. Velo, D. Zwanziger, Phys. Rev. 186, 1337 (1969). 



24 



[13] K. Johnson and E. C. G. Sudarshan, Ann. Phys. 13, 126 (1961). 

[14] W. K. Tung, Phys. Rev. Lett. 16, 763 (1966); Phys. Rev. 156, 1385 
(1967). 

[15] S. J. Chang, Phys. Rev. Lett. 17, 1024 (1966). 

[16] R. P. Feynman, M. Gell-Mann, Phys. Rev. 109, 193 (1958). 

[17] D. V. Ahluwalia, D. J. Ernst, Phys. Lett. B287, 18 (1992); Phys. Rev. 
C45, 3010 (1992); Int. J. Mod. Phys. E2, 397 (1993); D. V. Ahluwalia, 
T. Goldman, M. B. Jonson, Mod. Phys. Lett. A9, 439 (1994); D. V. 
Ahluwalia, Int. J. Mod. Phys. All, 1855 (1996); 

V. V. Dvoeglasov, Int. J. Theor. Phys. 37, 1915 (1998). 
[18] F. Sauter, Z. Physik 69, 742; 73, 547 (1931). 

[19] J. M. Jauch and F. Rohrlich, The Theory of Photons and Electrons 
(Addison- Wesley Publishing Co., Cambridge, Mass., 1955). 

[20] J. Schwinger. Phys. Rev. 82 664 (1951). 

[21] Ya. B. Zel'dovich, L. P. Pitaevskii, V. S. Popov and A. A. Starobinskii, 
Sov. Phys.-Usp. 14, 811 (1972) (Usp. Fiz. Nauk 105 , 781 (1971). 

[22] A. I. Nikishov, Sov. Phys.-JETP 30, 660 (1970) (Zh. Eksp. Teor. Fiz. 
57, 1210 (1969)); Nucl. Phys. B21, 346 (1970); Trudy FIAN 111, 152 
(1979). 

[23] M. S. Marinov and V. S. Popov, Sov. J. Nucl. Phys. 15, 702 (1972) 
(Yad. Fiz. 15, 1271 (1972)). 

[24] V. S. Vanyashin and M. V. Terent'ev, Sov. Phys.-JETP 21, 375 (1965) 
(Zh. Eksp. Teor. Fiz. 48, 565 (1965)). 

[25] W. Dittrich, J. Phys. A: Math. Gen. 11, 1191 (1978). 

[26] P. M. Lavrov, Izv. Vuz. Fiz. 3, 37 (1987). 

[27] A. Proca, Compt. Rend. 202, 1420 (1936). 



25 



[28] E. Petiau, Tesis, Paris, 1936. 

[29] R. J. Duffin, Phys. Rev. 54, 1114 (1938). 

[30] H. Kemmer, Proc. Roy. Soc. 173, 91 (1939). 

[31] I. M. Gel'fand, R. A. Minlos and Z. Ya. Shapiro, Representations of the 
Rotation and Lorentz Groups and Their Applications (Pergamon, New 
York, 1963). 

[32] M. A. Naimark, Linear Representations of the Lorentz Group (Perga- 
mon, New York, 1964). 

[33] Harish-Chandra, Phys. Rev. 71, 793 (1947). 

[34] W. J. Hurley and E. C. G. Sudarshan, Ann. Phys. (N. Y.) 85 , 546 
(1974). 

[35] F. I. Fedorov, Zh. Eksp. Teor. Fiz. 35, 495 (1958); The Lorentz Group 
(Nauka, Moscow, 1979). 

[36] N. Cabibbo, E. Ferrari, Nuovo Cim. 23, 1147 (1962). 

[37] A. A. Grib, S. G. Mamayev and V. M. Mostepanenko, Quantum Effects 
in Intensive External Fields (Moscow: Atomizdat, 1980). 

[38] S. I. Kruglov, J. Phys. G: Nucl. Phys. 21,1643 (1995); 22, 461 (1996). 

[39] W. Heisenberg and H. Euler, Z. Physik 98, 714 (1936). 

[40] V. Weisskopf, Kgl. Danske Videnskab, Selskab, Mat.-fys. Medd. 14, No 
6 (1936). 

[41] I. M. Ternov and V. G. Bagrov, Zh. Eksp. Teor. Fiz. 55, 2273 (1968). 

[42] V. N. Bayer, B. M. Katkov and V. M. Strachovenko, Yad. Fiz. 24, 379 
(1976). 



26 



